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Abstract. We investigate the quantum to classical transition of small inhomogeneous 
fluctuations in the early Universe using the decoherence functional of Gell-Mann and Har- 
tle. We study two types of coarse graining; one due to coarse graining the value of the 
scalar field and the other due to summing over an environment. We compare the results 
with a previous study using an environment and the off-diagonal rule proposed by Zurek. 
We show that the two methods give different results. 
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1) Introduction. 

The gravitational instability picture for galaxy formation assumes that the early Uni- 
verse started with a very smooth background on which small density fluctuations were 
superimposed. It is these small fluctuations which are ultimately responsible for the struc- 
ture in the present Universe. They have been amplified by the gravitational interaction 
since the beginning of the matter dominated era and produced the galaxies we see. 

In the sixties and seventies no theories were able to predict the existence of these 
perturbations, they were just postulated to be there. Zeldovich 1 and Harrisson 2 suggested 
that in order to fit the observation the initial spectrum of these perturbations must be 
roughly scale free. In 1980 Guth 3 proposed the Inflationary scenario to solve the horizon, 
flatness and monopole problems of the Big Bang. This scenario asserts that the Univserse 
went through a phase of very rapid expansion in its very early stage. The Universe would 
have expanded by a factor of at least 10 28 in a mere 10~ 32 seconds. 

It was soon realised that this very rapid expansion would have very interesting effect 
on fields especially the inhomogeneous part of the inflaton 4 ' 5 . The state of small inhomo- 
geneities undergoes parametric amplification during the inflationary period as soon as as a 
given mode crosses the Hubble radius. This lead to a scale free spectrum. It was therefore 
suggested that these inhomogeneities gave rise to the needed density fluctuations in the 
early Universe. 

It was argued that the quantum expectation value of the square of the field < (p 2 > 
can be interpreted as a statistical average of classical perturbations. The argument used 
by Guth & Pi 2 was that (< (p 2 >< ir"^ >) 1 ' 2 » % and thus quantum mechanical affects 
should be negligible. 

Interestingly enough this consideration is not invariant under linear canonical trans- 
formations. To see this more clearly the Wigner function can be calculated. In general 
the Wigner function is not positive and cannot represent a classical phase space density 
distribution but in the case of a gaussian state it is positive. So let's assume it can give us 
an idea of the classical phase space distribution. The Wigner function is defined as 




(1.1) 
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where p is the state of the system. Figure 1 depicts the 1 — a contour of the Wigner 
function for a mode k of a massless scalar field in the Bunch-Davies vacuum. Initially 
it is an ellipse rotating with frequency k/2n whose amplitude is adiabatic. As soon as 
the wavelength of the mode crosses the Hubble radius the ellipse stop rotating and gets 
elongated in the momentum direction. As seen from the picture the variance of cj) and tv^ 
are such that (< (p 2 >< >) x / 2 >> % but the surface of this ellipse remains %. Using a 

linear canonical transformation so that <p an d 7r</> are in the direction of the proper axis of 
the ellipse would give (< (fi 2 >< n 2 >) 1 / 2 = h. It is therefore difficult to understand why 
the quantum mechanical average can be substituted by a statistical one. 



Figure 1. 1 — a contour of a Wigner function for a mode of a scalar field which has 
crossed the Hubble radius. The area remains fi even if (< (ft 2 >< n 2 >) l t 2 >> fi. 

A lot of effort has recently been focused on understanding the transition between 
quantum and classical mechanics. It has been proposed that a measure of the classicality 
of a system is obtained by investigating the off-diagonal terms of the density matrix 6 ' 7 . 

In ref [8] such a criteria was used to investigate the classicality of the inhomogeneous 
quantum fluctuations in the inflationary period. It was shown that these fluctuations were 
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not classical if they were not interacting with an environment. A simple model of an 
environment represented by a single scalar field was constructed and it was shown that the 
off diagonal terms of the density matrix, in the configuration space basis, decreased rapidly 
as soon as the mode left the Hubble radius. The main problem with this approach is the 
assumption that when the off-diagonal terms in configuration space vanish, the system 
behaves classically. The density matrix gives information about the field at a given instant 
in time but it does not indicate how a small cell in phase space evolves. It tells us only how 
the sum of all these cells evolve. Classical behaviour requires each small cell of phase space 
to evolve independently of the others, that is, for there to be no quantum interference 
between different cells of phase space. 

In this paper we want to investigate a different approach to classicality, the one using 
the decoherence functional 9 ' 10 . After introducing the method in section 2 we investigate 
its consequences for perturbations in the inflationary universe. We study two types of 
coarse graining; one due to coarse graining of the value of the scalar field and the other 
by summing over an environment as in ref [8]. We compare the coherence length of the 
decoherence functional coarse grained from an environment with that of the density matrix 
and find striking differences. We discuss and conclude in section 3. 

2) The Decoherence Functional. 

In recent years an approach to the quantum to classical transition has considered not 
the state or eigenvalues of operators at a given time but rather focused on histories defined 
by a series of the value of fields at a different time 9 ' 10 ' 11 ' 12 . The idea is that a necessary 
condition for a system to be thought of as classical is that the probability sum rule for 
different histories should be obeyed. In other words interference should vanish. In such a 
case Griffith called them consistent histories. 

The tool to calculate the probability of an history is the decoherence functional. The 
fine grained decoherence functional for histories defined by the positions at all times can 
be defined through the path integral 

D[h\ h] = 8{q' f - q f ) exp i [S[q\ V )} - S[q( V )}} p(q' t , q t , Vl ) (2.1) 

where S is the action for the given history. In order for two histories h', h to be consistent 
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the decoherence functional must be diagonal. Except for very special cases, fine grained 
histories will not decohere. A possible way to get consistent histories is to coarse grain 
them. 

A coarse graining of this decoherence functional can be obtained by looking at histories 
with approximate position or momenta, or by summing over some field which is considered 
an environment. It is the latter case which corresponds to the the decoherence studied by 
Zeh 7 and Zurek 8 . A coarse graining can be defined as 

D c [h',h}= [ Vq' [ Vq6(q' f - q f ) e W i[S[q\ V )} - %fa)]]p(«U,»fc). (2.2) 
Jh> Jh 

The path integral over q(rf) is over all paths that start at qi at rji, pass through the intervals 
A 1 (?7i), A 2 (?7 2 )..., A n (?7 n ) at rji,r}2...r} n and wind up at qf at time rjf. Similarly for q'(rj) 
which goes through primed interval but end at the same endpoint qf. 

We will evaluate (2.2) for a scalar field evolving in the early universe both for coarse 
graining of the field or of an environment. A crucial question is how to model this envi- 
ronment. Any realistic model will be very complicated and hard to analyze. However, the 
basic physics should emerge from the simplest models. Hence we use a model which can be 
solved exactly: the system is a real massless scalar field <&i, (the inflaton), the environment 
is taken to be a second massless real scalar field $2 interacting with $1 by their gradients. 
This will permit us to compare our results with the ones in ref. [8] . We consider the fields 
in the de Sitter phase of an expanding Universe with scale factor a(t) = exp(Ht), where 
H is the Hubble constant. 

The action of system and environment is 

1 = J d^xy/g^dd^) 2 + (d^ 2 ) 2 + 2c(0 M $ 1 0"$ 2 )) (2.3) 

where g is the determinant of the background metric with line element given by 

ds 2 = a 2 {-drf + dx 2 ). (2.4) 

c is a constant measuring the strength of interaction between system and environment. 
We shall normalize the conformal time r\ such that r\ ranges between —00 and and 
a = —(Hi]) -1 with H~ l being the Hubble radius. 
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We will study two types of coarse graining. The first one will consist in summing 
over the field $2 which mimicks the environment. The second one will consist of coarse 
graining the value of the field $1. 

Our Lagrangian is quadratic in the derivatives of the fields and can hence be diago- 
nalized using fields <E> + and $_ for which the interaction term disappears. The coherences 
in the quantum state between $ + and are only given by the initial conditions. For 
example, we could choose an initial state where these coherences vanish. In this case, a 
pure state gives rise to a pure state reduced density matrix when summing over one of 
the fields. Decoherence of one field cannot occur by summing over the other one. We, 
however, suppose that the inflaton and the environment do not form the diagonal basis. 
This assumption is reasonable since any inflaton field (whose reduced density matrix we 
want) will interact with gravitational perturbations (part of the environment). 

We can expand the fields in harmonics in a box of fixed comoving volume (physical 
volume a 3 ) and investigate a particular wavenumber k = (k 2 + k 2 + k 2 ) 1 / 2 . As there is no 
coupling between modes with different k, we can consider a single wavelength and drop 
the index k for convenience. The Lagrangian reduces to 

L(q,q,r,r, V ) = ^-[q 2 + f 2 - k 2 q 2 - k 2 r 2 + 2c(qr - k 2 qr)]. (2.5) 

For simplicity we will consider histories described by only two values of <?'s, the value at 
time rji and rjf. When a hamiltonian exists we ca n rewrite the decoherence functional in 
the operator formalism as 

D{h u h 2 ) = Tr[U( Vf - Vi )PZ iriP (vi)PZ 2 r 2 UHv f ~ Vi)P^ trf ]- ( 2 - 6 ) 

In order to compare with the results using the density matrix, we will look for histories 
where r is considered as an environment and at first q is fine grained. The two histories 
that we will consider are the define by starting at either qi or q 2 at rji and ending up at 
qf. The decoherence functional becomes 

D(h u h 2 ) = N J dr 1 dr 2 dr f K*(qf,rf,r]f;q2,r2,r]i)K(qf,rf,r]f;q 1 ,r 1 ,r]i)p(q 1 ,r 1 ;q2,r2,r]i) 

(2.7) 
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and the propagator K is 



K (Qf,r f ,Vf,qi,ri,Vi) = 







2ttH 2 x_ 





(q 2 +r 2 + 2cq f r f )y f ( q 2 + r 2 + 2cq 1 r 1 )y. 



H 2 V 2 



+ 



H 2 Vl 



+ 



2k 3 (qfqi + TfT\ + cq/ri + cr/Qi) 
W 2 



where 

x = —k 2 i]frii sin kA + kA cos kA — sin kA 
yf = — k 3 i]fi]i cos kA — k 2 7]f sin/cA 
yi = —k 3 r]frii cos kA + k 2 r]i sin kA 
and A = rjf — t]i. If we assume that the initial to be 

•0(g, r, r]i) = a exp -b[q 2 + r 2 + 2aqr] 

where b is complex and a is real, we find that 



(2.8) 

(2.9a) 
(2.96) 
(2.9c) 



(2.10) 



D(h\, /12) = -Dexp 



where 



i(l - c 2 ) 
2x 



+ 



H 2 



exp - {q(A + q 2 A* + qi q 2 B) 

(2.11) 



A =[b 2 (l - a 2 ) + 66* (1 + c 2 - 2ca)]/(6 + 6*) 
B = - 266*(c-a) 2 /(6 + 6*) 



(2.12) 



D =aa* 



n 



1/2 fc 3 



6 + 6*/ 2ttH 2 x' 



This decoherence functional predicts a certain coherence length L^f, which is the maximum 
length (in configuration space squared) between histories over which interference is not 
exponentially suppressed. However to get a better measure of the decoherence of the 
decoherence functional, Djf, we should divide L^f by the probability width of the system 



Pdf. This is obtained by setting qi = q 2 in (2.11) and finding the length in configuration 
space squared where the probability not exponentially suppressed. We find that 

n - Pd f _ A + A*-B _ Abb*{c-a) 2 

~ L d f ~ A + A* + B ~ (i — a 2)(p + 6 *)2 • ^ A6 > 

When Ddf >> 1 we have significant histories decoherence. This measure for histories 
decoherence can be compared to the one used in ref [8] , using the off-diagonal terms of the 
reduced density matrix which is given by 

f ( n \ 1/2 

pred(qi,q2,Vi) = J i>(qi,r,rji)i>*(q2,r,rji)dr = aa* ( - - - j exp - (q 2 f + q\j* + qiq 2 g) 

(2.14) 

where 

g = -2a 2 bb*/(b + b*), f = (b 2 (l -a 2 ) + bb*) / (b + b*) . (2.15) 
In this case an analogous measure used was 

n ^ Pdm _ f + r-g _ (b + b) 2 -a 2 (b-b*) 2 , 91fi . 
dm ~ L dm ~ f + r + g- (l-a*)(b + b*)i ' { } 

This expression was analysed in ref[8] for the Bunch-Davies initial condition which corre- 
sponds to 

k 2 

a = c, b= ^— -. (2.17) 



The limit of interest is long after Hubble crossing where \krji\ « 1 which implies that 
b ~ k 3 /(2H 2 ) — ik 2 /(2H 2 rji). In this limit we see that D dm » 1, however D d f = 1. 
Thus we have a situation where an arbitrarily large decoherence of the configuration space 
density matrix corresponds to a maximally coherent decoherence functional. The Bunch- 
Davies vacuum is the ground state. If we perturb the initial coupling away from c then 
we will have some histories decoherence as well as density matrix decoherence. In this 
case the propagator will generate an imaginary contribution to a. This imaginary part 
will modify (2.13) and (2.16) in a way that cancels the divergence which would otherwise 
occur for a ^ c (a real) in the limit kt]i — > 0. In this case there may be closer relationship 



between the two decoherence measures. We can get an idea of the relative strengths of the 
two decoherence measures (for real a) by taking their ratios. We find 



L dm _ A + A*-B cjc - 2a)bb* 

L df f + f*-9 (re6) 2 + « 2 (zm6) 2 ' 1 ' } 

We can see that if c = the two coherence lengths agree. This can be easily seen form 
(2.7). The integral over r$ will be proportional to 8{r\ — T2) and thus the decoherence 
functional is proportional to the initial density matrix. It is rather surprising however that 
turning on the interaction from c = to c = 2a will increase the coherence of the decoher- 
ence functional relative to the density matrix. However for c > 2a the coherence length 
of the density matrix is larger than the coherence length of the decoherence functional. 
This shows that in this case there is no obvious correlation between the two decoherence 
measures and that the decoherence of the density matrix does not imply the decoherence 
of the decoherence functional or vice-versa. 

Another possibility to obtain decoherence is to coarse grain our system in configuration 
space, to which we turn now. In this case the histories are not defined by precise values 
of q but by a range determined by a (a variance) around a given value. The decoherence 
functional can then be obtained by integrating the fine grained one. The Ps are projectors 
on a range 2a of the fields. They are rather tedious to work with analytically. It will be 
useful to keep the analytical result simple so we use the gaussian pseudo-projectors 

P qiVi = j dzidnexp ^ ^ a i ~ ) r *)( r *» z *\ ( 2 - 19 ) 

They are not exactly projectors as 

P 2 ^ P (2.20) 

but are a sufficiently good approximation for our purpose. Substituting (2.19) into (2.6) 
we get 

D\h t M) = / dzidztdzsexp [- - - (Z3 -J f)2 ]D{ hl M) (2.21) 



where D(hi,h 2 ) is given by (2.11). We choose the Bunch-Davies initial condition (2.17) 
for (2.11). This is the most natural initial state to choose, it considerably simplifies the 
algebra and it ensures by virtue of (2.13) that any decoherence obtained will not be due 
to the environmental coarse grain. We can rewrite the result in terms of Q = q± + q 2 and 
5 = q± — q 2 and get 

D c (h u h 2 ) = iVexp[aiQ 2 + a 2 5 2 + a 3 qj + a 4 Q5 + a 5 q f 5 + a 6 Qq f ] (2.22) 

where 

-1 1 M + M* + 2V 
ttl ~2^2 + 4a 4 MM* -V 2 

-1 1 M + M*-2V 
a 2 =7r~o + 



a 3 = 



2a 2 4a 4 MM* - V 2 

V(M + M* - 2V) 
a 2 (MM* - V 2 ) 

M-M* 



(2.23) 



a4 



a 5 



o 6 



2(i 4 (M* - V 2 ) 

W X I 2 {M + M* -2V) 
a 3 (MM* - V 2 ) 

iV l / 2 {M - M*) 
a 3 (MM* - V 2 ) 



with 

1 2n,* ^(1-c 2 )^ (l-c 2 ) 2 k 6 a 2 

a 2 2xH 2 r}f 4x 2 H 4 

(2.24) 

_ (l-c 2 ) 2 fc 6 (7 2 
4x 

Investigating (2.22-24) shows that coarse grained histories that are determined by 
their approximate positions at various times are not exactly consistent. Exact decoher- 
ence is rather difficult to obtain so we investigate approximate decoherence. Histories are 
approximatively consistent if 



\ReD(h u h 2 )\ < eMin[ReD(h u h x ), ReD(h 2 , h 2 )\ 
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(2.25) 



This only means that the off-diagonal are much smaller than its corresponding diagonal 
part and thus the classical sum rules applies approximatively. e controls how good the 
approximation is. If we consider symmetrical histories (q± = —q2,Qf = 0) then it is easy 
to see from (2.22) that (2.25) translates mathematically as 



^« 1. 

a 2 



(2.26) 



We also want to be able to interpret the quantum mechanical average of operators as a 
statistical one. This implies that the coarse grain should be smaller than the fluctuations 
(Aq) 2 of the field. For the Bunch Davies vacuum (2.17), in the long after Hubble crossing 
limit, (Ag) 2 — > H 2 /k 3 hence we require 



a" << 



tf 2 



(2.27) 



In general (2.23) will be very long expressions. However they simplify greatly in the 
late time limit which implies that from (2.9a,2.9c) y$ — > —k 3 r]f and x — > ~ fc 3 3A (3^7// + A 2 ). 
We further consider the limit A — > 0, rji — > and rjf — > 0. We can take this limit while 
keeping an arbitrary constant proper time interval, St since ^ = —Hi]. In this limit we 



find that 



ai 



a 2 



a 3 



(24 



a 6 



-1 1 

+ j 



2a 2 
-1 



3_ {l-c 2 )^ 



if 2 



-l 



2a 2 

-(1-c 2 ) 
a 2 

a 5 -> 
2(1 -c 2 ) 



2if 2 + (1 - c 2 )fc 3 a 2 
3ii" 2 + (1 -c 2 )k 3 a< 



3_ + (l-c 2 )^ 



H 2 



-l 



(2.28) 



For our model in the late time limit (2.26) becomes 



(2.29) 



Equation (2.29) tells us that there is weak decoherence but not a significant amount. To get 
a better feel for this number it is worth comparing it to the long before Hubble crossing 
limit which gives a\/a2 = 1. Thus the after Hubble crossing limit does lead to some 
decoherence but not a significant amount. 

Assuming decoherent histories (5 = 0, q = Q/2) we find that (2.22) becomes, using 
(2.28) and (2.27) 



D c (h, h) « iVexp 



(2.30) 



Thus at late times the histories would be peaked about qf ~ q which is exactly the behavior 
of the classical motion. 

3) Discussion and conclusion. 

We see from (2.16) that the decoherence in the density matrix is due purely to the 
phase of the wave-function. This dependence on the phase is interesting since the phase 
can always be changed by a point transformation on the Lagrangian. We can see this as 
follows. A point transformation will transform the Lagrangian as 

L(q(t),m - L(q(t),m ~ J t f(Q(t),t) (3.1) 

which in turn means that the action transforms as 

S[q(t)] -> S[q(t)] - f(q f , t f ) + /(& U) (3.2) 

This point transformation doesn't affect the classical equation of motion because they are 
derived from the stationary action condition 8S = S[q(t)] — S[q(t) + 8q(t)] = where 
Sq(t) vanishes at the endpoints. However from the general expression U(qf,tf,qi,ti) = 
N Ylpaths e%S f° r t ne quantum propagator we can see that under the transformation (3.2) 
the quantum propagator transforms as 

U(q f ,t f ; qi, U) - e'^^U^ ,t f ; q i: U)e^^ (3.3) 

which in turn means that the wave function transforms as 

^(q,t) ^ e-^mt) (3.4) 
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Physics is generally considered invariant under the point transformation (3.1) because 
expectation values of functions of q and the physical momenta q are invariant (it is im- 
portant to remember that the canonical momenta does change with (3.2)). However the 
reduced density matrix of a subsystem in is not invariant to these point transformations. 
A point transformation is exactly what is being done when surface terms are dropped in 
a lagrangian. Using (3.3) and (3.2) we can see that the decoherence functional (2.7) is 
invariant under point transformations. This is an important difference between the two 
formalisms. 

In models with more general couplings we should expect decoherence of the reduced 
density matrix to depend not only on the phase but also on the real part of the exponent of 
the wave function. In this case there might be a simpler relation between the decoherence 
functional and the evolution of the density matrix. 

It is also interesting to investigate the influence functional for this model. Naively we 
might relate a diagonal decoherence functional with the existence of a noise kernel in the 
influence functional. Consider (2.2) where q — > (q, r),q' — > (q' ', r') and the r, r' coordinates 
are completely coarse-grained out. In this case (2.2) becomes 



D[h',h]= / Vq / Vq' / dq^dqfdq'f S(q f - q' f ) expz [S f [q( V )} - S f [q'( V )}} F[q( V ), q'(r,)} 
Jh> Jh J 

(3.5) 

where F [q(rj), q'{rf)] the influence functional is 



F[q(rj),q'(rj)} = J dndr^dr f dr' f 6(r f - r' f )p(q i ,r i ;q' i ,r' i ,r] i ) 
r( r fyf>Vf) 

/ VrVr'expi[S f [r( V )} + Si[q(r,),r(v)] ~ S f [r'(v)} ~ Si[q'(v), r'(v)}} ■ 

J (ri .r' .m) 



, , (3.6) 

( ( r f' r f^f) 

'(ri,r'i,Vi) 

For our model (2.5) with the Bunch-Davies initial condition (2.17) we find that the influence 
functional is 

~2 rVf ,>2 rvf 



F[q(rj),q'(ri)} = exp 



ic 
~2 



x exp 



(1 - c 2 )Mf - (1 - c 2 )b*q 



b f + b) 



(3.7) 
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A striking feature of (3.7) is the absence of a noise kernel that is typically associated with 
decoherence. This is due to the very special form of interaction we have chosen. The result 
is that in (3.5), there is no exponential suppression of widely separated histories and hence 
no histories decoherence. The influence functional will still not have a noise kernel even if 
the initial state does not have c = a. This shows that the absence of noise kernel does not 
imply a coherent evolution. 

We have shown that the decoherence functional shows some decoherence for the in- 
teraction given in eq.(2.3) for a wide selection of initial states. We have also shown that 
there is a surprising result for the case c = a as we have already mentioned. This case 
surely needs further study in order to understand why a mixed state can lead to a maxi- 
mally coherent (factorizable) decoherence functional. We also considered the possibility of 
decoherence after Hubble crossing though coarse graining the system field. We found this 
led to weak decoherence after Hubble crossing but probably not enough for an effective 
quantum to classical transition. 
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